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The collision of a high energy cosmic ray with a nucleon in the upper atmosphere could produce 
long-lived heavy particles. Such particles would be very penetrating, since the energy loss in matter 
scales as the inverse mass, and could reach a neutrino telescope like IceCube from large zenith angles. 
Here we study this possibility and focus on the long-lived stau of SUSY models with a gravitino LSP. 
The signal would be a pair of muon-like parallel tracks separated by 50 meters along the detector. 
We evaluate the background of muon pairs and show that any events from zenith angles above 80° 
could be explained by the production of these heavy particles by cosmic rays. 

PACS numbers: 13.85.Tp, 14.80.Ly 



I. INTRODUCTION 

The hierarchy problem has motivated an intense search 
for new physics during the past 20 years. Colliders like 
LEP, the Tevatron, or the B factories have explored the 
standard model (SM) at the quantum level but have not 
reached the energy or the sensitivity necessary to detect 
new physics. There is experimental evidence for neu- 
trino masses, whereas cosmological data strongly suggest 
the presence of a stable weakly-interacting massive par- 
ticle (WIMP) as the origin of the dark matter of the 
universe. However, these features could be easily added 
to the SM without changing its structure. Therefore, as 
we approach the search for extensions like supersymme- 
try (SUSY), technicolor or extra dimensions at the LHC, 
it is clear that we should never underestimate the SM. 

On the other hand, cosmic rays are another source of 
elementary particles of very high energy with the poten- 
tial to explore the physics beyond the SM. When a proton 
of 10 8 GeV from outer space hits an atmospheric nucleon 
it provides a center of mass energy y/s — \/2m^E around 
14 TeV. A small fraction of these protons (or of the sec- 
ondary particles with still enough energy) may then pro- 
duce exotic massive particles. Of course, the question is 
whether such an event could give any observable signal. 
In this paper we argue that this is the case, the long- 
lived charged particles present in some extensions of the 
SM could provide a distinct signature when they cross a 
neutrino telescope from large zenith angles. 

The process that we propose could take place in SUSY 
models with an exact i?-parity, a gravitino lightest SUSY 
particle (LSP) working as dark matter, and a charged 
next-to-LSP (NLSP) The 
collision of the cosmic proton with the nucleon could pro- 



duce any pair of SUSY particles, which would then decay 
promptly into the NLSP. Since the NLSP couples very 
weakly to the LSP gravitino, it will be long-lived and 
able to cross a kilometer-long detector like IceCube. The 
generic features of this framework could be also found in 
other extensions of the SM, like Little Higgs models [ill ]. 
These models may incorporate a T-parity [l2j separat- 
ing the standard and the exotic particles. The T-parity 
would forbid unobserved mixing between both sectors, 
tree-level four fermion operators, and would also make 
the lightest particle in the odd sector stable. If this parti- 
cle (constituting the dark matter) is very weakly coupled 
with the rest, the next-to-lightest one will be long-lived. 

The possibility to observe quasi-stable gluinos in 
IceCube [l]| has been considered in [3, HH in the 
framework of split-SUSY models with very heavy 
sfermions (liij . Here we will focus on non-colored par- 
ticles, which present some remarkable differences with 
the gluinos. In particular, as one of these particles propa- 
gates in matter it loses energy at a much smaller rate than 
an i?-hadron [13, EH ■ This fact makes it easier to confuse 
it with a muon, but it also lets the particle reach IceCube 
form larger zenith angles. To be definite we will consider 
a long-lived stau tr , although our arguments would be 
analogous for any massive charged particle: charginos, 
other sleptons, or vector like leptons that may appear in 
Little Higgs models. 

Other analyses of the production of exotic particle s by 
cosmic rays refer to primary neutrinos [13, H3, [HI, l23j |. 
Being weakly interacting, the relative effect of new 
physics on the neutrino-nucleon cross section may be 
larger (see below), however, it is difficult to make pre- 
cise predictions until the flux of cosmic neutrinos is de- 
termined. In contrast, our analysis here relies on a flux 
of primary cosmic rays that is well-known in the relevant 
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We will assume in the following that the t^j life-time is 
much larger than the propagation time through the Earth (see 
e.g. El). 
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FIG. 1: pN cross section to produce SUSY particles for differ- 
ent values of the stau and the (common) squark/gluino mass 
and a proton energy of 10 8 GeV. 



FIG. 2: Cross section to produce any pair of colored SUSY 
particles in a pN and a ttN collision for different incident 
energies and a squark/gluino mass = 150, 300 GeV. 



energy region 10 5 -10 8 GeV. 

The outline of the paper is as follows. In Sec. HH we 
calculate the rate of long-lived stau pairs produced by 
collisions of primary and secondary cosmic rays in the 
atmosphere. We discuss then in Sec. IIIII the signal of 
these pairs at IceCube and the background from muon 
pairs. Section [IVI includes a summary of our results. 

II. PRODUCTION OF STAU PAIRS 

At energies above 10 4 GeV the decay length of nucle- 
ons, charged pions and kaons is much larger than their 
interaction length in the air. Therefore, as they propa- 
gate in the atmosphere the probability that one of these 
hadrons (h) collides with a nucleon (N) to produce new 
physics is just Q 

A rr hN 

-P h x(E)~^. (1) 

In this expression a 1 ^ 4 is the cross section to produce 
the exotic particle(s) X and <j!^ a the total cross sec- 
tion of the hadron with the air (we assume A = 14.6 
nucleons in a nucleus of air and neglect nuclear ef- 
fects). The (default) cross sections with the air used by 
CORSIKA IH above 10 4 GeV can be approximated by 
4 a « Cj + Cf ln(E;/GeV) + C£ ln 2 (£7GeV), with the 
constants given in Table U Since a!^ a is above 100 mb, it 
is apparent that this probability will be very small and 
that it would be much larger for a neutrino propagating 
in matter. 

The cross section to produce SUSY particles in a 
hadronic collision depends basically on their mass and 



on whether they have strong interactions. All the cross 
sections at the parton level can be found in Ref. [26| . 
Collider bounds on SUSY particles with prompt decay 
into a neutral LSP are around 250 GeV for gluinos 
and squarks, and 100 GeV for the stop, the sbottom, 
charginos, and charged sleptons These bounds, how- 
ever, may not apply if the particles decay instead into a 
long-lived charged or colored SUSY particle. For exam- 
ple, in order to minimize the SM background a recent 
analysis [28| of jets plus missing momentum at the run II 
of the Tevatron imposes a veto on events with an isolated 
electron or muon with large transverse momentum. How- 
ever, gluino or squark events will include here final staus 
(instead of neutralinos) that could be taken by isolated 
muons. We are not aware of specific bounds on the col- 
ored SUSY spectrum in stau NLSP models. Notice that 
bounds based on the delay in the time of flight versus 
a muon or the anomalous ionization of the staus should 
also be specific, as they are based on the absence of slow- 
moving (/3 < 0.6) charged particles, but here the staus 
get an extra boost from the decay of the parent squark or 
gluino. Through the paper we will then consider slepton, 



TABLE I: Constants defining the total cross section with the 
air o% a « C# + C? ln(£/GeV) + C h z ln 2 (S/GeV). 



h 


CH [mb] 


C? [mb] 


Cl [mb] 


N 


185.7 


13.3 


0.08 




100.5 


16.9 


0.00 


K 


79.7 


13.9 


0.05 



3 



"' — I 



10" 



I I 1 1 ■ ■ ■ 1 1 — I I I III III — I I I ■ I ■ ■ 1 1 — I I 1 1 1 1 1 1 1 — I I 1 1 llllj — I I 1 1 1 IN 1 



mx = 150 GeV 
300 GeV 




I I I I I ■ I ■ ■ 



10 4 10 5 10 6 10 7 10 8 10 9 10 10 
E [GeV] 



I 



10 
1 b- 

10" 1 L 

10" 2 

10" 3 t- 
10 -4 L 



10" 



I I I I Mill — 'III I Mill 1 I I I Mill 1 I I I Mill 1 — I I I 1 1 hi 



300 GeV - 1 




' '"I " 



10 4 10 5 10 6 10 7 10 8 10 9 
£W [GeV] 



FIG. 3: Average angle between the two staus (in the labora- 
tory frame) for = 150, 300 GeV. 



FIG. 4: Flux of stau pairs produced by cosmic rays in terms 
of the (total) stau energy E^t on production for m% = 150, 
300 GeV and r\ = 1. We include the flux of muon pairs 
produced with an opening angle above 10~ 4 rad. 



chargino and neutralino masses as low as 100 GeV and 
colored SUSY particles above 150 GeV. 

In Fig. Q] and [2] we plot the total hN (h = p, n) 
cross sections to produce pairs of these SUSY particles 
for different SUSY masses (left panel) and for values 
of the hadron energy between 10 4 and 10 11 GeV (right 
panel). We have used the CTEQ6M H (MRSS O) 
parton distribution functions for baryonic (mesonic) in- 
teractions, with the renormalization scale /i — 0.2m^ 
suggested by a NLO calculation [3l| . We observe that 
the cross section to produce directly a pair of long-lived 
staus of 100 GeV is much smaller than via the produc- 
tion and prompt decay of colored particles of mass around 
200 GeV. In the latter case, we estimate that the final 
stau will carry a fraction n 

m% + m% ~ 
^Hfe- (2) 

of the energy of the parent gluino or squark. The average 
angle of the stau pair (in the lab frame) for different 
energies of an incident proton are given in Fig. [3] 

To evaluate the production rate of stau pairs by cosmic 
rays we need the total flux of hadrons: primary plus sec- 
ondary nucleons, pions and kaons produced at any depth 
in the atmosphere and with enough energy to create staus 
in the collision with an air nucleon. This analysis has 
been carried out in djj assuming a flux of primary nu- 
cleons 

d$w _ / E \ ~ 2 ' 7 nucleons 
dE ~ ■ \1 GeV J cm 2 s sr GeV ' 1 ' 



for energies up to 10 6 GeV, with a spectral index that 
changes to 3 in the interval 10 6 -10 9 GeV and goes back 
to 2.7 at higher energies {e.g. [32[). We will use in the 
following the results on the flux of secondary hadrons 
(nucleons, pions, and kaons) derived there. 

The flux of quasi-stable staus produced via the prompt 
decay of a pair XX' of SUSY particles is then 

E f ^V^,{E). (4) 

h=N,n,K Jli ™™ 

In Fig. d] we plot the differential flux (d$/d£Vr) of stau 
pairs produced by cosmic rays for squarks and gluinos 
masses of 150 and 300 GeV and r\ = 1. 



III. BACKGROUND OF MUON PAIRS AND 
SIGNAL AT ICECUBE 

The flux of stau pairs produced high in the atmosphere 
needs to propagate down to the core of IceCube, about 
two kilometers under the antarctic ice, to be observed. 
In addition, the possible signal faces a strong background 
of muon pairs crossing the detector. We plot in Fig. 0] 
the flux (d^/d-E^) of muon pairs 2 produced by cosmic 
rays of energy Eh > 10 4 GeV. We include only the events 
where the two muons are produced with an opening angle 



2 We neglect the muons from tau decays as they are a 1% cor- 
rection to this flux. 
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FIG. 5: Range of staus and muons of different energy. 
The dashed lines show the integrated column depth of the 
Earth from the center of IceCube for increasing zenith angles 
(Aftnax = 10°). 
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FIG. 6: Integrated number of stau and muon pairs with a 
minimal separation of 50 meters at the detector for different 
values of a minimal zenith angle. 



above 10 rad, since smaller angles imply a separation 
between the two muon tracks that can not be resolved 
at IceCube (see below). Notice that this requirement 
cuts off muon pairs with an invariant mass near thresh- 
old, Vs < 1 GeV, where the PDFs are mostly unknown 
and the process would be better described in terms of 
hadronic resonances. 

The propagation of muons and heavy charged particles 
in matter is well understood. For a muon of energy > 
2to m the mean energy loss per column density (measured 
in g/cm 2 ) can be approximated as 

dE 

- = «m + / W - (5) 

where a M s» 2 x 10~ 3 GeV cm 2 /g describes ioniza- 
tion effects and ^ « 4 x 10~ 6 cm 2 /g accounts for 
bremsstrahlung, pair production and photohadronic pro- 
cesses. The solution of Eq. ([5} provides an approximation 
of the total range of (initially) very relativistic muons, 
which we consider in the following. 3 

For a stau, at the lowest order ionization effects co- 
incide (a^ ~ a^) whereas the other effects depend 
mainly on the velocity of the particle, which implies 
/3? w /3„m„/my. For sw 100 GeV, this means that 
a stau of energy above 10 5 GeV losses 10 3 times less en- 
ergy than a muon of the same energy, but below 500 GeV 
they deposit energy at a similar rate. In our analysis we 



3 Below E as 2m the ionization energy loss grows like 1//3 2 (see 
e.g. [27t). 



will use the approximation for the range of a stau pro- 
vided in [T3| and will neglect losses through electroweak 
interactions (l8j . as they are not important for the stau 
energies that we obtain. 

In Fig. [5] we plot the range of staus and muons of en- 
ergy between 10 3 and 10 9 GeV. We give the correspon- 
dence between integrated column depth of the Earth (see 
e.g. [331 ]) and zenith angle for several trajectories ending 
at the center of IceCube. We observe, for example, that 
whereas a muon of E = 10 7 GeV has a range of about 
25 km water equivalent (w.e.) and can reach IceCube 
from a zenith angle max ~ 86°, the range of a stau of 
the same energy is 360 times larger, which makes it able 
to reach IceCube from zenith angles of up to 6* max w 105°. 

Another relevant observable is the separation of the 
two particles when they cross IceCube. This depends on 
their angle at the creation point (see Fig. [3]) and the dis- 
tance from that point to the telescope. The interaction 
length of a 10 7 GeV proton in air is around 4 g/cm 2 (its 
cross section is a^ a w 400 mb), which corresponds to 
an altitude of about 20 km in the atmosphere. There- 
fore, if a primary proton creates a stau pair it will do it 
around that altitude. The production of a stau pair by 
a secondary hadron will typically occur along the second 
interaction length, finishing at around 15 km, and so on. 
To estimate the distance between the parallel stau tracks 
at IceCube we will assume that they are created at a 
height H w 15 km. This implies that stau pairs coming 
from zenith angles of 60° , 80° , and 100° fly an approxi- 
mate distance of about 30, 90, and 2300 km, respectively 
to reach the center of IceCube. 

In Fig. [5] we compare the number of stau and muon 
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pairs reaching IceCube with a separation larger than 50 
meters, so that the two tracks can be resolved [34|]. We 
plot the flux of these particles coming from zenith angles 
larger than the value given in the x-axis (e.g. 0min = 
corresponds to pairs coming from any direction). Wc 
show the cases — 150, 300 GeV and r] — 1, 0.7. 
We observe that from zenith angles between 80° and 95° 
there is a possible signal with no background from muon 
pairs. 

IV. SUMMARY AND DISCUSSION 

Cosmic rays may be continuously producing massive 
particles when they collide with nucleons in the upper 
atmosphere. If these particles are long-lived they will be 
able to reach a detector like IceCube, about two kilome- 
ters under the ice. To have a sizeable production rate (or- 
der 1 per year and square kilometer) the particles should 
be produced through strong interactions. This would be 
the case for a quasi-stable stau resulting from the prompt 
decay of a gluino or a squark. We have studied in some 
detail the possibility to observe such an event. 

The heavy staus would be produced in pairs at al- 
titudes around 15 km, and as they approach IceCube 
the two staus would separate. In principle, they could 
be confused with a muon pair: a muon and a stau of 
500 GcV would give in IceCube a very similar signa- 
ture. We have shown, however, that it is possible to 
reduce the dimuon background below the signal. Above 
E ss 500 GeV muons lose energy in ice much faster than 
the staus. As a consequence, while muons will never 
reach IceCube from directions close to the horizon, staus 
can come from zenith angles of up to 110°. In addition, 
larger zenith angles mean also larger distance of flight 
and, in turn, larger separation between the two tracks at 
IceCube. 



We obtain that any events with two muon-like tracks 
separated by more than 50 meters coming from zenith an- 
gles of 80°-100° would be a clear signal of heavy charged 
particles produced by cosmic rays high in the atmosphere 
(see Fig. 16]). 

One may wonder if this signal could also be distin- 
guished from possible stau pairs produced by primary 
neutrinos, which have been extensively considered in the 
literature. We find two main differences with such an 
event. Suppose that the primary comes horizontally 
(9 = 90°). In our case the staus will be created high 
in the atmosphere, whereas in a neutrino event the in- 
teraction to produce them will typically take place deep 
inside the ice. This implies a smaller distance between 
the two stau tracks along IceCube. A second important 
difference is that while neutrino events could come basi- 
cally from any direction, the staus produced by primary 
protons vanish at zenith angles above about 115°. 

We think that the study of inclined (9 > 60°) two- 
muon events at IceCube would be interesting by itself. 
For example, it can be used as an indirect measure of the 
total (primary plus secondary) flux of hadrons of 10 4 - 
10 9 GeV and, thus, as a test for the different codes that 
simulate extensive air showers. Going to zenith angles 
below the horizon the SM background vanishes, leaving 
some room for exotic physics. 



Acknowledgments 

We would like to thank Francisco del Aguila, Fernando 
Cornet, Andreas Ringwald, and Christian Spiering for 
useful discussions. This work has been supported by 
MEC of Spain (FPA2006-05294) and Junta de Andalucia 
(FQM-101). We also acknowledge financial support from 
a grant CICYT-INFN (07-10). 



[1] A. Brandenburg, L. Covi, K. Hamaguchi, L. Roszkowski, 
and F. D. Steffen, Phys. Lett. B617, 99 (2005), hep- 
ph/0501287. 

[2] W. Buchmiiller, K. Hamaguchi, M. Ratz, and 
T. Yanagida, Phys. Lett. B588, 90 (2004), hep- 
ph/0402179. 

[3] K. Hamaguchi, Y. Kuno, T. Nakaya, and M. M. Nojiri, 

Phys. Rev. D70, 115007 (2004), hep-ph/0409248. 
[4] J. L. Feng and B. T. Smith, Phys. Rev. D71, 015004 

(2005) , hep-ph/0409278. 

[5] K. Hamaguchi and A. Ibarra, JHEP 02, 028 (2005), hep- 
ph/0412229. 

[6] O. Cakir, I. T. Cakir, J. R. Ellis, and Z. Kirca, (2007), 

hep-ph/0703121. 
[7] R. H. Cyburt, J. Ellis, B. D. Fields, K. A. Olive, and V. C. 

Spanos, JCAP 0611, 014 (2006), astro-ph/0608562. 
[8] J. R. Ellis, A. R. Raklev, and O. K. Oye, JHEP 10, 061 

(2006) , hep-ph/0607261. 

[9] J. R. Ellis, K. A. Olive, Y. Santoso, and V. C. Spanos, 



Phys. Lett. B588, 7 (2004), hep-ph/0312262. 

[10] J. L. Feng, S. Su, and F. Takayama, Phys. Rev. D70, 
075019 (2004), hep-ph/0404231. 

[11] N. Arkani-Hamed, A. G. Cohen, E. Katz, and A. E. Nel- 
son, JHEP 07, 034 (2002), hep-ph/0206021. 

[12] H.-C. Cheng and I. Low, JHEP 09, 051 (2003), hep- 
ph/0308199. 

[13] IceCube, J. Ahrens et al, Astropart. Phys. 20, 507 

(2004) , astro-ph/0305196, |http : //icecube .wise . edu/| 
[14] J. I. Illana, M. Masip, and D. Meloni, Phys. Rev. D75, 

055002 (2007), hep-ph/0611036. 
[15] J. L. Hewett, B. Lillie, M. Masip, and T. G. Rizzo, JHEP 

09, 070 (2004), hep-ph/0408248. 
[16] N. Arkani-Hamed and S. Dimopoulos, JHEP 06, 073 

(2005) , hep-th/0405159. 

[17] M. H. Reno, I. Sarcevic, and S. Su, Astropart. Phys. 24, 

107 (2005), hep-ph/0503030. 
[18] Y. Huang, M. H. Reno, I. Sarcevic, and J. Uscinski, Phys. 

Rev. D74, 115009 (2006), hep-ph/0607216. 



6 



[19] G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 

(1999), hep-ph/9801271. 
[20] I. Albuquerque, G. Burdman, and Z. Chacko, Phys. Rev. 

Lett. 92, 221802 (2004), hep-ph/0312197. 
[21] M. Ahlers, J. Kersten, and A. Ringwald, JCAP 0607, 

005 (2006), hep-ph/0604188. 
[22] X.-J. Bi, J.-X. Wang, C. Zhang, and X. Zhang, Phys. 

Rev. D70, 123512 (2004), hep-ph/0404263. 
[23] I. F. M. Albuquerque, G. Burdman, and Z. Chacko, 

(2006), hep-ph/0605120. 
[24] J. I. Illana, M. Masip, and D. Meloni, (2006), hep- 

ph/0612305. 

[25] D. Heck, G. Schatz, T. Thouw, J. Knapp, and 
J. N. Capdevielle, CORSIKA: A Monte Carlo code 
to simulate extensive air showers, Report FZKA- 
6019 (1998), Forschungszentrum Karlsruhe; http://www- 
ik.fzk.de/corsika , FZKA-6019. 

[26] S. Dawson, E. Eichten, and C. Quigg, Phys. Rev. D31, 



1581 (1985). 

[27] Particle Data Group, W. M. Yao et al, J. Phys. G33, 1 
(2006). 

[28] DO, V. M. Abazov et al, Phys. Lett. B638, 119 (2006), 

hep-ex/0604029. 
[29] J. Pumplin et al, JHEP 07, 012 (2002), hep-ph/0201195. 
[30] P. J. Sutton, A. D. Martin, R. G. Roberts, and W. J. 

Stirling, Phys. Rev. D45, 2349 (1992). 
[31] W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas, 

Nucl. Phys. B492, 51 (1997), hep-ph/9610490. 
[32] Particle Data Group, S. Eidelman et al, Phys. Lett. 

B592, 1 (2004). 
[33] R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, As- 

tropart. Phys. 5, 81 (1996), hep-ph/9512364. 
[34] M. Ribordy, Nucl. Instrum. Meth. A574, 137 (2007), 

astro-ph/0611604. 



